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We present a fully-exclusive next-to-next-to-leading order (NNLO) QCD calculation for vector 
boson production in hadron-hadron collisions. The calculation is implemented in a parton level 
Monte Carlo program, which includes y — Z interference, finite-width effects, the leptonic decay 
of the vector bosons and the corresponding spin correlations. The code allows the user to apply 
arbitrary (though infrared safe) kinematical cuts on the final-states and to compute distributions 
in the form of bin histograms. We show some illustrative numerical results at the Tevatron and 
the LHC. 
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Vector boson production in hadron collisions, the well known Drell-Yan (DY) process [[]]], has 
a special role for physics studies at hadron colliders. Having large production rates with relatively 
simple experimental signatures, this process is important for detectors calibration, it gives strin- 
gent informations on Partem Distribution Functions (PDFs), is a strong test for perturbative QCD 
predictions and may signals effects from physics beyond the Standard Model. 

It is therefore essential to have accurate theoretical predictions for the vector-boson production 
cross sections and related distributions, which are predicted by perturbative QCD as an expansion in 
the strong coupling as- The next-to-next-to-leading order (NNLO) QCD corrections (i.e. £?(a|)) 
have been calculated analytically for the total cross section [§] and the rapidity distribution of the 
vector boson The fully exclusive NNLO calculation has also been performed [Q, g]. Further- 
more, electroweak corrections up to 0(a) have been computed for both W [§|] and Z production |J7|]. 

The computation of higher-order QCD corrections to hard-scattering processes is a hard task. 
Difficulties arise from the presence of infrared singularities at intermediate stages of the calculation 
that prevent a straightforward implementation of numerical techniques. For the above reason, fully 
exclusive cross-sections in hadron collisions have been computed so far only for Higgs boson 
production [§, ^, 10, 11] and the Drell-Yan process [|], 

In this contribution we present a recent fully exclusive NNLO QCD calculation for vector 
boson production in hadron collisions The calculation uses the NNLO extension of subtraction 



formalism introduced in Ref. [10]. The method is valid in general for the production of colourless 



high-mass systems in hadron collisions. 

We consider the hard-scattering process: 



h l +h 2 ^V(q)+X, 



(1) 



where the colliding hadrons hi and h 2 produce the vector boson V (V = Z/y*,W + or W ), with 
four-momentum q and high invariant mass \[q*-, plus an inclusive final state X. 



Following Ref. [JlOj, we observe that, at LO, the transverse momentum qj of V is exactly zero. 
This means that, as long as qj ^ 0, the (N)NLO contributions are given by the (N)LO contributions 
to the final state V + jet(s) [O]: 



aa (N)NLO\<lT^0 



ao (N)LO 



(2) 



We compute do^o^ by using the subtraction method at NLO [[D], [l4|] and we treat the remain- 
ing NNLO singularities at qj = by the additional subtraction of an universal 1 counter-term 
d®(N)w constmcte d by exploiting the universality of the logarithmically-enhanced contributions 
to the transverse momentum distribution 2 . Schematically we have 



(N)NLO 



y {N)NLO 



do, 



^.y+jets 



(N)LO 



do, 



(N)W 



(3) 



where ^Z\ NLO is a process-dependent coefficient function necessary to reproduce the correct nor- 
malization 



'it depends only on the flavour of the initial-state partons involved in the LO partonic subprocess. 
2 For the explicit form of the counter-term see Refs.[|lo[ [l5|]. 
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We have encoded our NNLO computation in a parton level Monte Carlo event generator. The 
calculation includes finite-width effects, the y — Z interference, the leptonic decay of the vector 
bosons and the corresponding spin correlations. Our numerical code is particularly suitable for the 
computation of distributions in the form of bin histograms. 

In the following we present some illustrative numerical results for Z and W production at the 
Tevatron and the LHC. We consider u,d,s,c,b quarks in the initial state. In the case of pro- 
duction, we use the (unitarity constrained) CKM matrix elements V lic i = 0.97419, V us = 0.2257, 
Vub = 0.00359, V cd = 0.2256, V cs = 0.97334, V cb = 0.0415 from the PDG 2008 @. We use 
the so called scheme for the electroweak couplings, with the following input parameters: 
Gf = 1.16637 x 10~ 5 GeV 2 , m z = 91.1876 GeV, T z = 2.4952 GeV, m w = 80.398 GeV and 



TV = 2.141 GeV. As for the PDFs, we use MSTW2008 [JTSj] as default set, evaluating a s at each 
corresponding order (i.e., we use (n + l)-loop as at N"LO, with n = 0, 1,2). We fix the renormal- 
ization (Hr) and factorization (jip ) scales to the mass of the vector boson my. 
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Figure 1: Rapidity distribution of an on- shell Z boson at the LHC. Results obtained with the MSTW2008 
set (left panel) are compared with those obtained with the MRST2004 set ( right panel). 



We start the presentation of our results by considering the inclusive production of e + e~ pairs 
from the decay of an on-shell Z boson at the LHC. In the left panel of Fig. [I] we show the ra- 
pidity distribution of the e + e~ pair at LO, NLO and NNLO, computed by using the MSTW2008 
PDFs[18]. The corresponding cross sections are Olo = 1.761 ±0.001 nb, Gnlo = 2.030 ±0.001 nb 
and Gnnlo = 2.089 ±0.003 nb. The total cross section is increased by about 3% in going from NLO 
to NNLO. In the right panel of Fig. [T] we show the results obtained by using the MRST2002 LO 



3 



NNLO QCD correction to vector boson production at hadron colliders 



Giancarlo Ferrera 



[19] and MRST2004 [20] sets of partem distribution functions. The corresponding cross sections 
are a w = 1.629 ±0.001 nb, o NL0 = 1.992 ±0.001 nb and a NNL0 = 1.954 ±0.003 nb. In this case 
the total cross section is decreased by about 2% in going from NLO to NNLO. 
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pp W + +X -> Yv+H 



MSTW2008 NNLO 
JR09VF NNLO 
ABKM09 NNLO 




Figure 2: Rapidity distribution of an on-shell W + boson at the Tevatron. The NNLO result obtained with 
the MSTW2008 set are compared with those obtained with the JR09VF and ABKM09 sets. 

We next consider the production of an on-shell W + boson at the Tevatron. In Fig. ^| we show 
the rapidity distribution of the W + at NNLO, computed by using the MSTW2008 PDFs [JTs|] . We 
also show, for comparison, the NNLO prediction by using the JR09VF [|TJ and ABKM09 @ 



PDFs. The corresponding total cross sections are Gnnlo^ = 1-349 ±0.002 nb, o^^[p = 1.338 ± 
0.002 nb and o^™^ = 1.391 ±0.002 nb. The differences between the three results can be reach, 
in the central rapidity region, the level of about 5%. 

We finally consider the production of a charged lepton plus missing pj through the decay of 
a W boson (W = W + ,W~) at the Tevatron. The charged lepton is selected to have pr > 20 GeV 
and | tj | < 2 and the missing pj of the event is required to be larger than 25 GeV. The transverse 



mass of the event is defined as mj = y2/?' r p™ lss (l — cos0), where is the angle between the 
the pj of the lepton and the missing pj. In Fig. |5| we show the transverse mass distribution at 
LO, NLO and NNLO: the accepted cross sections are Olo = 1.161 ±0.001 nb, Onlo = 1.550 ± 
0.001 nb and Gnnlo = 1.586 ±0.002 nb. Since at LO the W boson is produced with zero transverse 
momentum, the requirement pj lss > 25 GeV sets my > 50 GeV. As a consequence at LO the 
transverse mass distribution has a kinematical boundary at mj = 50 GeV. Around this boundary 



there are perturbative instabilities due to (integrable) logarithmic singularities [23]. We also note 
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Figure 3: Transverse mass distribution for W production at the Tevatron. 

that, below the boundary, the NNLO corrections to the NLO result are large. This is not unexpected, 
since in this region of transverse masses, the ff((Xs) result corresponds to the calculation at the first 
perturbative order and, therefore, our ^(cuf) result is actually only a calculation at the NLO level 
of perturbative accuracy. 

We have presented a fully exclusive NNLO QCD calculation for vector boson production in 
hadron-hadron collisions. Our calculation is directly implemented in a partem level event gener- 
ator. This feature makes it particularly suitable for practical applications to the computation of 
distributions in the form of bin histograms. For illustrative purpose, we have shown some selected 
numerical distributions at the Tevatron and the LHC. 
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